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Abstract 



We study supersymmetric structure of the seh-isospectral crystahine chains formed 
psj ■ by copies of the mutuahy displaced one-gap Lame systems. It is generated by the 

N[N — 1) integrals of motion which are the first order matrix differential operators, 
by the same number of the nontrivial second order integrals, and by the N third order 
Lax integrals. We show that the structure admits distinct alternatives for a grading 
operator, and in dependence on its choice one of the third order matrix integrals plays 
r"! ! either the role of the bosonic central charge or the role of the fermionic supercharge 

to be a square root of the spectral polynomial. Yet another peculiarity is that the 
■ set of all the second order integrals of motion generates a nonlinear sub-superalgebra. 

^ . We also investigate the associated self-isospectral soliton chains, and discuss possible 

^ \ physical applications of the unusual extended supersymmetry. 

^" 

o. 

1 Introduction 

Quantum periodic finite-gap systems are associated with completely integrable models. Any 
such a system is characterized by a proper nontrivial integral of motion that is a higher order 
Lax operator P [1]. In the infinite period limit they transform into refiectionless (soliton) 
systems. Their band structure is effectively encrypted in a Burchnall-Chaundy operator 
identity [21 E] that can be presented in a form of a supersymmetric-like relation [3] 

{V,V} = 2P{H), (1.1) 

where P{H) is a spectral polynomial and if is a Hamiltoniar0. 



^Because of a higher order nature of the Lax integral, the indicated relation formally corresponds to a 
higher-derivative generalization of the supersymmetric quantum mechanics [5] , whose construction naturally 
arises after truncation of the parasupersymmetric quantum mechanics [6]. Nonlinear supersymmetry of 
an extended system with higher order superscharges can be related to the Crum-Darboux transformations 
[3 IHl ini HD] in the way like the usual supersymmetric quantum mechanics is related to the first order Darboux 
transformations In the unextended case, however, relation (|l.ip has a nature of a hidden nonlinear 

supersymmetry |12j . in which the Z2-grading is provided by a refiection operator |13j that identifies the 
parity-odd Lax integral as a fermionic supercharge [4] . In the extended case of the chains of finite-gap systems 
we study here, the Crum-Darboux and hidden supersymmetric structures naturally meet. For the earlier 
treatment of finite-gap systems in supersymmetric and related contexts see also [TH [TSl [HI HTJ [T51 [TH [5D] . 
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Because of the presence of the nontrivial integral V, a super symmetric extension of 
finite-gap systems has a more rich structure in comparison with that of a usual N = 2 
supersymmetric quantum mechanics. An example of a physical system with such an unusual 
supersymmetry is provided by a model of a non-relativistic electron in periodic magnetic 
and electric fields |21j. In a generic case of the N = 2 superextension of an n-gap system by 
means of a Crum-Darboux transformation of the order m, 1 < m < 2n + 1, where 2n + 1 
is a number of band edge states, supersymmetric finite-gap system is characterized by two 
further supercharges of the order {2n + 1 — m) [21] . The anti-commutator of supercharges 
of the orders m and (2n -\- 1 — m) generates the Lax integral of the order 2n + 1. Another 
important peculiarity is that such superextended systems admit various choices for the grad- 
ing operator. For some of them, matrix diagonal Lax integral takes a role of one of the 
supercharges that annihilates all the 2(2?t, + 1) band edge states of the extended system. 
Unlike the anti-diagonal supercharges of the orders m and {2n -|- 1 — m), it also distinguishes 
the left- and right- moving Bloch modes inside the allowed bands |22]- This is similar to a 
role played by a momentum operator for a free particle that can be treated as a zero-gap 
system. 

Sometimes, a supersymmetric partner happens to be just a spatially displaced initial 
finite-gap periodic or non-periodic system. The Witten index takes then a zero value for 
such a self-iso spectral system [T7] even in the case of unbroken supersymmetry [1 ^ [T^ 123] . 

Recently, it was observed [22] that a phase transition between the kink-antikink and 
kink crystalline phases in the Gross-Neveu model [211 |25l |26l [271 l2H] is accompanied by the 
structural changes in the associated supersymmetric self-isospectral one-gap periodic Lame 
system. Any of the two first order supercharges of the latter can be taken as the Bogoliubov- 
de Gennes Hamiltonian in the Andreev approximation, in which superpotential is identified 
with a gap function (a condensate field). The first order supercharges are constructed from 
the Darboux displacement generators of the associated second order Schrodinger system. Yet 
another peculiarity of such a system is that the first order Bogoliubov-de Gennes Hamiltonian 
possesses its own, exotic hidden nonlinear supersymmetry [22]. A certain infinite period limit 
applied to the one-gap self-isospectral system reproduces either the supersymmetric structure 
of the Dashen, Hasslacher, and Neveu kink-antikink baryons [29] as a Darboux dressed form 
of a free massive Dirac particle [30] , or superextended version of the Callan-Coleman-Gross- 
Zee kink solution [291 EI] of the Gross-Neveu model [22] . 

A pair of the second order supercharges of the N = 2 superextended one-gap Lame 
system is generated by a sequence of the two first order Darboux displacements, while the 
Lax operators of the mutually shifted Lame subsystems are produced by closed third order 
Crum-Darboux loops. The second order supercharges and the third order Lax integrals 
include a dependence on an auxiliary, virtual displacement parameter. It is natural to try 
to extend the model by taking N > 2 mutually displaced copies of Lame systems. This 
produces then the question: 

• What supersymmetric structure will be generated in such an extended Lame system 
by considering the higher order unclosed (m > 2) and closed (m > 3) sequences of 
Darboux displacements? 

In the present paper we study the supersymmetric structure of the self-isospectral crystalline 
chain of mutually displaced one-gap Lame systems, and investigate in the same context the 
associated self-isospectral non-periodic (soliton) chains of reflectionless Poschl- Teller systems. 
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The paper is organized as follows. In the next section we construct the first order Dar- 
boux displacement generators for the one-gap Lame system. In section 3 we apply them to 
construct self-isospectral chains of Lame systems, and investigate the higher order general- 
izations of the Darboux displacement generators. In section 4 we study a general structure of 
the exotic nonlinear supersymmetry of the self-isospectral crystalline iV-term chain. Section 
5 is devoted to the discussion of alternative choices for the Z2-grading operator admitted by 
such an extended supersymmetric structure. In section 6 general theory is illustrated by an 
example of the = 3 crystalline chain. In section 7 we consider the infinite period limit 
that produces a self-isospectral non-periodic chain and its supersymmetric structure. In the 
last section we conclude with a discussion of possible physical applications of the revealed 
unusual nonlinear supersymmetry. 



2 Darboux displacement generators: periodic case 

Consider a one- dimensional self-adjoint Schrodinger Hamiltonian H{x) = + with a 
periodic potential V{x). Require that the system admits a family of the first order Darboux 
displacement generators V{x; A) = ^ + (f{x; A), 

V{x;X)H{x) = H{x + X)V{x;X), (2.1) 

which depends on a continuous parameter A. Then it can be shown that V{x) has to be a 
one-gap Lame potential [HI |19]. The Jacobi form of the one-gap Lame system is [32] 

H(x) = -— + 2k^sn'^x - k\ (2.2) 
dx^ 

A modular parameter k, < k < 1, fixes a real, 2K, and imaginary, 2iK', periods of the 
potential. Here and in what follows we do not indicate explicitly the dependence of elliptic 
and related functions on k, and use a notation i = for imaginary unit. The chosen 

value of the additive constant fixes the level of the lower edge of the valence band to be 
zero, and the one-gap spectrum of (12. 2p is <j{H) = [0, /c'^] U [l,oo), where < fc' < 1 is 
a complementary modular parameter, k'"^ = 1 — k"^. The infinite period limit (A; — )■ 1 ^ 
2K — )■ oo, 2iK' — )■ ivr, snx — )■ tanhx) of (12. 2p corresponds to a reflectionless Poschl- Teller 
system with one bound state in the spectrum, while in another limit A; — )■ 0, (12. 2p reduces to 
a free particle system. 

To construct a one-parametric Darboux displacement generator, we discuss shortly some 
properties of Lame system (12. 2p . Solutions of the stationary equation if\I'(a;) = E'^lx) are 
given by the Bloch functions of the form 

(a;) = ^^Q^^^ exp [TxZ{a)] , E = dn^a, (2.3) 

where H, 9 and Z are the Eta, Theta and Zeta Jacobi functions [32], [33]. Under translation 
for the period, they transform as 

^^(x + 2K) =exp[Ti2K/s:(a)]*^(x), /€(«) = - iZ(a) , (2.4) 
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where K{a) is a quasi-momentum. Energy E is given here as a function E(a) = dn^a of 
a complex parameter a. This is an elhptic function with the same modular parameter k, 
and its period parallelogram in complex plane a G C is a rectangle with vertices in 0, 2K, 
2K + 2iK', and 2iK'. On the border of the indicated period parallelogram, function dn a 
takes real or pure imaginary values, and so, E is real. The vertical sides a = i/3 + K, 
< /3 < K', and a = if3, < (3 < K', correspond, respectively, to the valence, < E < k''^, 
and the conduction, 1 < E < oo, bands, where the quasi-momentum K{a) is real. The 
horizontal sides a = iK' + /3 and a = /3 with < /3 < K correspond to the prohibited bands 
—oo < E < and k''^ < E < 1, where K{a) takes complex values. Inside the allowed bands, 
(12. 3p are the two Bloch modes propagating to the left (the upper index) and to the right (the 
lower index). On the edges of the bands, they reduce to the standing waves described by 
the periodic, dnx = dn (a; + 2K) {E = 0), and antiperiodic, cnx = — cn [x + 2K) {E = k'"^) 
and sna; = — sn(a; + 2K) {E = 1), functions. 

Like the ground state dnx (« = K + iK', E = 0), non-physical Hamiltonian eigenstates 
in the lower prohibited band — oo < < are nodeless functions, which can be used 
to construct a one-parametric family of the first order Darboux generators. As dn {—u) = 
dn (m + 2K) = — dn ('u + 2iK') = dnu, it is convenient to introduce a notation a = — 2r + iK', 
and assume that r G M while keeping in mind that E — oo for r — )■ riK, n G Z. By 
shifting the argument, x ^ x + r, for the wave function (12.31) with the upper index we get 
^;2^+'^'(a; + r) = c{t)F{x; r), where c(r) is a nonzero x-independent multiplier, and 

|||^exp[x3(r)]=F(x;r). (2.5) 

Here 

3(r) = Z(2. + iK') + i^ = Z(2.) + Hi|^^ ,2.6, 

is, up to a factor — i, a quasi-momentum of the Bloch state (12. 5p . 3(r) = — iK(— 2r + iK'), 
that is an odd function of r. A nodeless function F{x; r) is quasi-periodic in x, F{x + 
2K; r) = exp[2K3(r)]F(a;; r), periodic in r, F{x;t + 2K) = F{x]t), and for x 7^ it 
undergoes infinite jumps from to -|-oo at r = riK, n E 1^. It also satisfies the relations 
F{x;—t) = F{—x;t) = 1/F(x;r). In the case of the ground state 3(— K/2) = 0, and 
function (12. 5p reduces to a periodic function (/c')~^/^dn(x — y). 
Let us consider now a first order differential operator 

V{x; r) = F{x- r)-^— ^ = - A(x; r), V^{x; r) = -V{x; -r) , (2.7) 
ax r [x; T) ax 

whose zero mode is F{x;t), V{x;t)F{x;t) = 0. Function A(x;r) = F'{x;t)/F{x;t), 
F'{x]t) = ^F{x;t), reads 

A(x;r) = 3(r) + Z(x-r)-Z(x + r) 

cn 2rdn 2r , 2 / \ / \ /r^ o\ 

= \- k sn 2rsn(x — r )sn(x + r) . (2.o) 

sn2r 

It obeys the Riccati equations 

A2(x; r) ± A'(x; r) = 2khn^{x ±T)-k'^ + £(r) , (2.9) 



4 



where 

e(r) = -E{-2t + iK') = cn22r/sn22r . (2.10) 
Another important property is that the following three-term linear combination, 

A(x; r) + A(x + r + A; A) + A(x + A; -r - A) 

= + 3(A) + 3(-T - A) ^ A) , (2.11) 

is x-independent. The function g{T, A) possesses the symmetry properties g{T, A) = g{X, r) = 
g{T, —A — r) = —g{—T, —A) and can be presented in a form 

, 1 - cn2rcn2Acn2(r + A) 

gi'T, A = ^ . ■ 2.12 

' ^ sn2rsn2Asn2(r + A) ^ ' 

From Eqs. ([2lD and ^J}, we find a relation A(^; r) = ^(^(^ - r - iK'); r), i.e. by (12TT|) 
A(^;r) can be presented by a three-term sum of quasi- momenta taken at three values of 
arguments which sum up to the zero. We also will need the identity 

A'{x + T + \;\) - A(x + A; r + A)A(x + r + A; A) + g{T, A)A(x; r) = 
-^(A^(a;;r) + A'(x;r) + 5(r)), where 6{t) = 1 + k"^ - 3sn-^2T , (2.13) 

which follows from (ES]) and f l2TT]l . 

By the Riccati equations (12.91) . the operators (12. 7p factorize the Lame Hamiltonian (12. 2p . 

P'^(x;r)P(a;;r) =i7(a; + r)+£(r), ©(x; r)r'^(s; r) = if(x - r) + e(r) , (2.14) 

with (I2.10p playing a role of a factorization constant. 

Note here that the second eigenstate vl/2^^+'^ (x + r) of the shifted Lame Hamiltonian 
operator (12. 2p reduces, up to an inessential r-dependent multiplier, to F{x + 2r; — r), that 
is coherent with a relation E{~2t + iK') = E{2t + iK'). In correspondence with this 
observation, a change r — )■ — r in the first relation from (I2.14p and a subsequent shift x — )■ 
a; -|- 2r transform this first factorization into an equivalent form V{x + 2r; r)P^(x + 2r; r) = 
H{x + r) + £:(t), that is just the second relation from (I2.14p with the argument x shifted for 
2r. 

From (12.140 it follows that (12. 7p are the sought for Darboux displacement generators, 

V{x;t)H{x + t) = H{x-t)V{x;t), V\x]t)H{x - t) = H{x + t)V\x]t). (2.15) 

As V\x] t) = —V{x; — t), it is sufficient to consider only the first intertwining relation from 
(I2.15P while the second follows from it via a simple change r — )■ — r. 



3 Chains of one-gap Lame systems 

In this section we construct higher order generalizations of Darboux generators, that will 
lead us naturally to the chains of Darboux-displaced one-gap Lame systems. 
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A mutual spatial Darboux displacement between the two systems in fl2.15p is 2r while 
their 'average coordinate' is x. For generalization it is convenient to characterize each of the 
two related systems by its own shift parameter by introducing the notations 

1 1 

Tab = ^{n - To) = -Ua , Xab = X + -{Ta + u) = Xba- (3.1) 

Then Xab + Tab = x + Tb, Xab — Tab = X + Ta, and relations ( I2.14p . (I2.15P can be presented in 
the form 

'DabVl = -VabVba = Ha + Sab, (3.2) 

T^abHb = HaVab, (3.3) 

where we have introduced further notations 

Vab = V{Xab;Tab) = -Vl^, Ha = H{x + Ta), S ab = £{Tab) = Sba , (3.4) 

see Fig. 1. Since the superpotential A, the Darboux displacement generator Vab and the 
factorization constant Sab blow up at Tab = "^K, n G Z, we suppose that Tab 7^ '^K. 




Figure 1: a) The Darboux displacement generator Vab transforms eigenstates of Hb into 
those of the translated system Ha, see ( 13. 3p . b) The Hamiltonian Ha as a sequence of the 
two Darboux displacements (13. 2p . 

Making use of the relation ( 13. 3 p , one can define the second order operator 

Bab/X = VaxVl = -VaxVxb , B^^^ = Bba/X , (3.5) 

where as for Tab, we assume that Tax, T\b 7^ nK. Like the first order operator Dab, it inter- 
twines the same two systems Ha and Hb, 

Bab/xHb = HaBab/X, (3.6) 

via a chain of the two Darboux displacements, Bab/xHb = —VaxT^xbHb = —VaxH\D\b = 
—HaVaxl^xb = HaBab/X- In this chaiu, there appears an intermediate system H\, which from 
the viewpoint of our pair of basic systems Ha and Hb is of a virtual, auxiliary nature. To 
stress a virtual nature of the displacement parameter A, we indicate it in a special way (with 
slash) in notation for the second order Crum-Darboux intertwining operator B. From ( 13. 5 p 
we find a relation 

Baa/X=VaxVl^ = Ha+EaX- (3.7) 

Hence, Bab/x is a kind of a non-Hermitian generalization of the Lame Hamiltonian operator. 
In correspondence with (13.71) . the second order intertwiner Bab/x, unhke Vab, is well defined 
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(for Tax, T\f) 7^ riK) also in the case when Tab = nK.. The virtual parameter on the right-hand 
side in (13. 7p appears only in the additive term. We also have 



^ab/X^lh/x = ^ab/X^ba/X = {Ha + eax){Ha + ^bx) 



(3. 



Making use of (12 .13^ . we find that a specific linear combination of the second, i3, and the 
first order, D, intertwining operators, 

yab = —^ah/X — dabxT^ab , 3^]b = 3^6a 5 (3.9) 

does not depend on the virtual parameter A, where we have introduced a notation 

gabX = girab, TXa) = d^Tab) + ^{Ux) + d{TXa) (3.10) 

for a function of displacement parameters, which is completely antisymmetric in the indices, 
dabx = —Qbax = —QaXb- Note the cyclic order of the indices on the r.h.s. of (I3.10p . The 
explicit form of the intertwining operator yab = y{xab; Tab), yahHb = Hayab, is given by 



(P A( ■ ) 



k SB. {x + t) + sn ^2r . 



(3.11) 



y{x;T) 
From (13. 9p we get also 

Sab/X = }3ab//i + {gabfi — gabx)T^ab , (3.12) 

that corresponds to a change of the virtual displacement parameter, see Fig. 2. The co- 
efficient in fl3.12p before the first order intertwining operator has a cyclic representation in 
terms of the quasi-momentum, Qacfji - QacX = 3(^aA) + ^(^Ac) + di^cf,) + 3(^Ma), cf. (|3.10p. 





+ (9abv~9abx) 



a 




Figure 2: The second order intertwiner as a sequence of the two Darboux displacements; 
second line corresponds to a change A — )■ of the virtual parameter f l3.12p . 



Intertwining operators of the first and the second order allow us to construct a nontrivial 
integral for Lame system Ha, 



'Pa — T^abyba + ^ab^ab — yab'P'ba — ^abC. 



ab ■ 



(3.13) 



[Va,Ha] = 0, where Va = V{x + Ta), Cab = C^Tab) = -Cba, C{t) = (^(r, ^K) = 
dn2r/(sn2rcn2r). Integral fl3.13p is nothing else as the Lax operator for one-gap Lame 
system ( 12. 2p . whose exphcit form is 

p(x) = — + (1 + - 3fc2sn\)— - 3A;^sna;cnxdnx, = -V . (3.14) 

\AjtAj KajJLj 
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Relation ( I3.13p can be presented in the equivalent form 
where 

^abc — ^abi^dabc ^ab^ • 

Making use of the equivalent representation for the function (12. lip . (12.120 . 

C(r)5(r) - C(/i)e(/i) 



(3.15) 
(3.16) 



£(r) - e{fi) 



one can check that the three-index object defined in (13.161) possesses the same antisymmetry 
properties as gabc, ^abc = -^bac = -^acb- We can also write 



Pa — —^ab/J^ba + QabcHa + iabc — 

see Fig. 3. With the help of (I3.17P and relations 

^ab^bc^ca S, 



-P'ab^ba/c 



QabcHa ^abc i 



(3.17) 



abc 



^ab^ac ^~ ^ab^bc ~l~ ^ac^bc "^dabc^abc 



£ab + £ac + ^bc — 9. 



abc 



0, 



(3.18) 



we find that the Lax integral and the Hamiltonian satisfy the Burchnall-Chaundy operator 
identity 

-V^ = P{H), P{H) = H{H-k''^){H-l), (3.19) 

where P{H) is a spectral polynomial of the one-gap Lame system (12. 2p . One can show that 
in correspondence with (I3.19p . the physical states (12.30 are also the eigenstates of the Lax 
operator, V'^'^{x) = ^ir|{E{a))^y P{E{a)), where ri{E) = —1 for the valence and -|-1 for 
the conduction bands [22]. The Lax integral distinguishes the left- (^E'") and the right- (\E'°) 
moving Bloch modes inside these bands, and annihilates the band-edge states. 




Figure 3: Two representations for the third order (Lax) integral: Pa 

\Bav/\P'ua = —\P'av'Pv\P'\a — \P'a\P\vP'va- 



' 2 P^au^ua/X 



Proceeding from the definition of the second order operator Bab/\ as a composition of the 
first order Darboux displacement generators, one can generalize the picture by treating the 



8 



intermediate system if a not just as a virtual one but on equal grounds as the systems Ha 
and Hb. 

The second order intertwining generator can be generalized for the case of the third order. 
Making use of the already known lower order relations and identity fl3.12p . one finds 

T^abDbc^cd = —^ac/bT^cd = —T^ab^bd/c 

= —{Ha + €dc)T^ad + {gacd — gacb)Bad/c = —{Ha + £ab)T^ad + {Qbda — gbdc)Bad/b ■ (3.20) 

The variation in the form of the two last expressions is in the intermediate indices. Such 
an ambiguity (freedom) related to the intermediate indices also appears in the higher order 
generalizations of the intertwining relations and integrals (Darboux loops). If we put d = a 
in (13. 20 p . the triple product of the Vs reduces to a third order Darboux loop (I3.15p . while 
in the two last expressions the first order operator V and the coefficient before the second 
order intertwining operator are singular. A careful treatment of these expressions in the 
limit sense d — )■ a reproduces correctly the r.h.s. of the relation (I3.15p . 
In the next, the fourth order case, we get in the same way 

T^abT^bc^cdT^df = ^ac/b^cf/d = {Ha + Edc + Qadc — gadf)^af/d 

+ {gacd - gacb) {Ha + S fd)Vaf. (3.21) 

Taking in the last relation f = a, we find that the fourth order integral for Ha can be 
presented as follows 

T^abT^bcT^cdT^da = {Ha + edc){Ha + Sad) + {gacb — gacd){Va + gacdHa + ^acd) , (3.22) 

i.e. it reduces to a function of Ha and Va- Continuing, we find that in the general case the 
closed (loop) sequence of the n Darboux transformations Vab^'Db^b^ . . .T>b„_^a is an integral 
of motion for the system Ha of the form 

Vab,Vb,b, . . . Vb„_,a = h{Ha) + h2{Ha)Va , (3.23) 

where hi 2{Ha) are certain polynomials of Ha- An unclosed sequence of n Darboux transfor- 
mations reduces, analogously, to a combination of the first and the second order intertwining 
operators with coefficients to be some functions of the intertwined Hamiltonians Ha and He-, 

Vab,Vb,b, . . . Vb„_,c = h{Ha)Vac + f2{Ha)Bac/b„_, = VaJl{H,) + J2 (^^c) • (3.24) 

Index bn-i in Bac/b„-i can be changed for any other intermediate index by employing identity 

(EISD. 

We conclude therefore that the higher order open chains of the Darboux displacement 
transformations reduce as differential operators to linear combinations of the two basic 
blocks: the first, V, and the second, B, order Darboux displacement generators with co- 
efficients to be certain functions of the intertwined Hamiltonians. In the case of the closed 
(loop) chains, they reduce, analogously, to a linear function of the third order Lax integral 
V with coefficients depending on the Hamiltonian. No new structures do appear in addition 
to these sets of the first and the second order intertwining generators and the third order 
Lax integrals, which will play a role of the basic blocks in the associated supersymmetric 
construction, to the discussion of which we pass over in the next section. 
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4 Supersymmetry of self-isospectral periodic chains 



In this section we introduce a kind of the A^-extended system to be a self-isospectral chain 
of one-gap Lame systems, and study general characteristics of the supersymmetric structure 
associated with it. 

Consider a chain of > 2 one-gap Lame systems which we describe by a matrix Hamil- 
tonian 

n = diag{Hi,...,HN). (4.1) 

Here we use the same notations as in (13.41) . Ha = H{x + Ta), a = 1, . . . , N, and assume that 
the set of the shift parameters is restricted by the condition Tab = ^{Tb — Ta) 7^ nK for 
any pair of indices a b, i.e. we suppose that the arguments of Hamiltonians of any two 
subsystems are shifted mutually for any distance to be different from the real period n2K. 

Introduce a symbol iab defined by Cat = +1 for a < 6, = — 1 for a > b, and iab = if 
a = b. We imply that (cab)^""^^ = 4?,, n e Z, while (eab)^'' = +1 for a 7^ 6 and (cafe)^" = 
for a = b. We also introduce the N x N matrices: 

(Sf ),^. = {eabfiStS'; + S',S]) , (Sf )^^. = ieabiStS'; - S'.S]) , (4.2) 
i^t),, = eab{5t5] - 5\5]) , (r^)^^. = {eabfm^, + 5\5]) . (4.3) 

ab „„J irab 



In (14.21) and (14. 3 p we assume that a ^ b, and so, the first factor in definition of S" and f 

^b 



can be omitted. By definition, all the four matrices are symmetric in indices a, b, S^^ = S^"^ 



k = 1,2, 3, I" = I . For N = 2 they reduce to the three Pauli and the unit matrices, and 



for > 2 satisfy the same algebra Sf Sf = 5ijF^ + ieijk^f, i,], k = 1,2, 3. 

Making use of the intertwining relations from the previous section, we find that the 
system (14.11) is characterized by the A^(A^ — 1) nontrivial integrals 

{Sf),^ = (i)'+' (Sf )^^. {e.,yv,, , 1 = 1,2, (4.4) 

which are the matrix differential operators of the first order, and by the same number of the 
integrals of the second order 

{Qt/x),^ = {^15) + 5\5^^){ie.,r'B,,/^ , 1 = 1,2, (4.5) 

[n,Sf] = Q, [n,Q^l]=Q. (4.6) 

In definition of the integrals ( 14. 5p we suppose that the virtual parameter t\ can take inde- 
pendent values for each pair of indices a ^ b, and for any of the two values of the index /; 
the only restriction, as before, is Ta\,T\b 7^ nK. On the other hand, relation (I3.12p means 
that any second order integral Q with the changed value of the virtual parameter is a linear 
combination of the initial operator Q and of the first order integral S. 

In accordance with the introduced notations, the integrals with the indices I = 1 and 
/ = 2 are related by 

St = iLfSt, Qti^ = iSf g?^, . (4.7) 

This is coherent with the fact that I"'' and YI^ are also the integrals of motion for "H, 
['H,!"''] = [7^,23^] = 0, and they act, respectively, as the identity and (J3 = diag{l,—V) 
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matrices in the two-term subsystem specified by the indices a ^ h. Also, the following 
relations are valid: 

j^ab ^ j^ba^ M'^^W"^ = when all a, b, c, d are distinct , (4.8) 

where any of M and W is 5*; or Qi/x. 
Matrix 

r = dtagil, -1, . . . , i-ir) , = 1 , (4.9) 

is an (zero order) integral of motion for the system (14. ip . [H, T] = 0, and can be taken as a 
grading operator. It identifies the integrals 5*"^ and Qf^x with a — 6 = 2?t, + 1 as fermionic 
operators, {F, M"-^} = 0, and those with a — 6 = 2?7, as bosonic, [T,M''^] = 0. To identify the 
superalgebraic structure generated by the integrals of motion, we compute anti-commutators 
between fermionic integrals (supercharges), and we take commutators between bosonic, and 
between bosonic and fermionic integrals. In accordance with (14.81) . corresponding commuta- 
tors, [M'^^,M'"^], and anti-commutators, {M'''', M"'^}, take zero values when all the indices 



a, b, c, d are distinct. Nontrivial anti-commutators for fermionic supercharges are 

{sf, s'^} = 2Si,nr\n + Sab) , (4.10) 

{Sf, S'::} = {eabYiet^r {SU-^YQT/b + (l - SiJcacQT/b) , (4.11) 

{Qt/x, Qt/x} = 25imr\n + Saxm + ebx) , (4.12) 

{Q^/x.Qt/x} = {eaby'\e,,r-' {Sim{-^y-'QT/x + (1 - Sim)eacQT/x) (n + e,,) , (4.13) 

{Qt/x, S:^} = 2ieaby~\e,arr' {5U-^)\9abxn + Ux) - (1 - 5vm)L) , (4.14) 

{Qt/x. S'^} = {eaby-\hcr {^U-iy^^eacST + (l - 5im)ST) in + etc) ■ (4.15) 
The nontrivial commutators are given by 

[Sf, S'^] = -2iieaby^\haril - 6lm)^fin + Sab) , (4.16) 

[5f , St:] = Keabyihcr {Slmi-iy^'eacQT/b - (1 - Slm)QT/b) , (4.17) 

[Qt/x, Qt/x] = -2i(e,,)'(e,„)— 1(1 - SiJ^fCH + SaxW + 6,x) , (4.18) 

[Qt/x, Qwa] = i{eaby-\e,cr^' {Sim{-^yeacQT/x + (1 - Sim)QT/x) (n + e,x) , (4.19) 

[gr;„ S^J;] = 2i(e„,)'(e,,)"^Sf {i-iy+'6i^L + (1 - 6,^)igabxn + Ux)) , (4.20) 

[Qt/x, S';:] = Keaby-\e,cr {SiU-^ysr + (l - 5i,n)eacSr) in + e^J . (4.21) 



In correspondence with the comment on the change of the virtual index we made above, 
without loss of generality we put the same value for it in both second order integrals in each 
of the corresponding (anti)-commutators in fl4.12p . (14.130 . fl4.18p and (14.190 . 
In (14.141) and (I4.20p there appears a nontrivial matrix operator 

L = -idiag{Vi,V2r-- ,Vn), (4.22) 
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composed from the third order Lax operators Va = V{x + Ta) = —^{T>abBba/x + Bab/xl^ba), see 
Eq. fl3.13p . The bosonic integral L is a central element of the superalgebra with the grading 
operator F, 

[L, n] = [r, L] = 0, [L, Sf] = [L, Q^/a] = . (4.23) 

For the self-isospectral chain system (14. ip . we have N obvious (bosonic for f l4.9p ) third order 
integrals = —idiag{0, . . . , 0, Va, 0, . . . , 0), the sum of which corresponds to the central 
element L. To write the commutation relations for them in the general case of the A^-terms 
chain, it is convenient to define the linear combinations of L*^, Lf' = I^^'L and Lg'' = ^f'L, 
remembering that for N > 2 not all them are linearly independent. All these third order 
integrals commute between themselves, while their nontrivial commutators with the first and 
the second order integrals are 





{eabr+'5i2i^t 




(4.24) 


I T ab obc] •m, 
[-^i 5 •^ml — 1 ' 






(4.25) 


r T ab /^ab ] 
i-'^l 5 Vm/Aj - 




™M/A, 


(4.26) 


[Lf, Q^/a] = -i' 


''^'{hat'eU 




(4.27) 



where / = 1,2, and we denote 

[LSr' = Qf,^{U + Sab) - eabSt{9abxn + Ux) , (4.28) 

[LQ]% = eabQt/xigabx'H + Ux) + St{H{H + 2k'^e-J - e^x - Sbx) + SaxEbx) • (4.29) 

Though on the right hand side of (14.281) . there appears explicitly a virtual parameter A, the 
complete combination of the integrals there does not depend on A. This can be checked by 
making use of definitions (14.41) . (14. 5 p and Eqs. (13. 9p and (13.160 . 

We see that the set of the first, the second and the third order integrals of motion, which 
are Hermitian matrix operators, generate a kind of nonlinear superalgebra. A nonlinearity 
is related to the fact that some of the (anti)-commutators of these integrals are quadratic in 
the Hamiltonian "H, or include "H or 'H? as a multiplier at other integrals. The results on the 
general form of intertwining relations from the previous section show that no new independent 
integrals do appear in addition to those we already found. The interesting property of this 
supersymmetric structure is also that the set of the second order integrals Ql^^ taken with 
the same value for the virtual parameter A (i.e. with the same shift parameter tx) together 
with the Hamiltonian l-L form a closed nonlinear sub-superalgebra, see Eqs. (14.120 . 14.13p . 
(I4.18P and (14.190 . The second order integrals can be reduced to a 'standard' form with a 
prescribed, (any) fixed value of the virtual parameter by means of relation (I3.12p . 

5 Alternative choices for the Z2-grading operator 

The choice (14. 9 p for the grading operator is not unique. The permutation of diagonal elements 
in (14. 9p . or multiplication of some of them by —1 changes the identification of integrals as 
fermionic and bosonic ones, and, as a consequence, some commutation relations will be 
changed for anti-commutation relations and vice versa. This does not change, however, the 
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bosonic nature of the diagonal matrix integrals , and a global conclusion on a nonlinear 
nature of superalgebra. 

There are alternative choices for the grading operator which involve reflections in the 
coordinate x and in the shift parameters. They provide some new features for the super- 
algebraic structure. Let us discuss some of such alternatives. Consider the reflection in x 
(parity) operator 7?., TZx = —xTZ, IZti = TiTZ, 7^^ = 1, and the operator T that reflects any 
of the shift parameters, including the virtual ones, Tr,; = —TiT, Tx = xT, = 1, and 
commutes with TZ. The product of these two operators is a nontrivial, nonlocal integral of 
motion for our chain system, [}ZT,7{\ = 0, and its square equals 1. So, it can be identified 
as another sort of the grading operator, 

T = nr. (5.1) 

The operator f anti-commutes with all the first, Sf^, and the third, Lf^, order integrals, and 
commutes with the second order integrals Qf/x- ^^^^^ '^^^^ "^""^ and L"^ are the fermionic 
operators, while are the bosonic ones. To compute the super algebraic structure for such 
a choice of the grading operator, we have to use coherently with the described identification 
of bosonic and fermionic generators the corresponding (anti)-commutators from the previous 
section, which should be supplied with the nontrivial anti-commutation relations that involve 



the third order integrals, 

{Lf, S:^} = 2i"^{ea,r+'6n{^fr[LSr' , (5.2) 

{Lt^S'^} = nhaY-'i^lTlLSr , (5.3) 

{Lt\Q^J:^,} = ~2i"^5n{^fr'''[LQ]%, (5.4) 

{^f ,Qwa} = -i'''^\e,ay-\^lT^'[LQ]%, (5.5) 

{Lf,L:^} = 2{J:fy^"^P{'H), (5.6) 

{Lf = (F" - eab^f){e,ay-\e,,r-'P{n) . (5.7) 



Under such a choice of the grading operator, the spectral polynomial of the chain, -P('H) = 
HCH — k''^){'H — 1), appears exphcitly in the superalgebraic structure. This choice, therefore, 
is coherent with the structure of a hidden, bosonized supersymmetry fll.ip [12] that is present 
in each of the chain subsystems Ha The relation (15.61) particularly reveals the property 
of the Lax operators that is essential for physical applications: each third order differential 
operator Li here is an annihilator of the three band edge states in the spectrum of the 
corresponding chain Lame subsystem Hi, see refs. [2T1 [221 ES] for the further details. In the 
case of the choice of the grading operator discussed in the previous section, this peculiarity 
of the Lax operators does not show up in the superalgebraic relations. 

The product of the two operators, (14.91) and (15.11) . can also be chosen as a grading 
operator. Yet other possibilities are associated with the introduction of the operators of 
the permutations of the displacement parameters, T"-^ = T^°-, defined by T'^^rf, = TgT"'^, 
T°-^Tc = TcT"-^, T'^^x = xT°-^. Combining such operators with the matrix structures S^^'gs, 
and reflection operators TZ and T, one can construct more integrals which can be taken as 
the grading operators. This does not add something essentially new to the structures we 
already observed, and we do not discuss these other possibilities here. 
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In the last section we present some further arguments in favor of necessity to consider al- 
ternative choices for the grading operator (alongside with the choice discussed in the previous 
section) in the context of possible physical applications. 



6 Supersymmetric structure of the = 3 chain 



In the case of the two-term chain {N = 2), we have a, 6 = 1, 2. A complete set of independent 
integrals is formed by the two integrals of the first order, SP, I = 1,2, the two integrals of the 
second order, QlJ^, and by the two integrals of the third order, = L and L^^. In the list 
of the (anti)-commutation relations of the integrals there do not appear (anti)-commutators 
which involve the generators with three different indices a, b, c. For the discussion of the case 

= 2 we refer to [22]. Here we consider in more detail the next case = 3 to illustrate 
explicit matrix form of the involved structures. 

The A^ = 3 Hamiltonian is 

fH, \ 

n = I H2 . (6.1) 

The system has six trivial matrix integrals of motion, 

= I 1 I , I^^ = I I , I^^ = I 1 I , (6.2) 
\0 0/ 







^0 







)' 





1 


• 




^0 








= -1 , S^^ = , Sf = 1 , (6.3) 



which appear explicitly in the superalgebraic (anti)-commutation relations. This set contains 
only three linearly independent matrices. Three nontrivial integrals to be the first order 
differential operators are 





f 




sr = 1 









{ 







= I I , 5f = I -P23 I , (6.4) 
and the other three are S'2'' = iTif'S^^. Six second order integrals of motion are given by 







Ql% = 1321/x , Q\% 



.B 









^13/A 











'31/A 











^0 





° \ 










^23/A 






^32/A 


/ 



(6.5) 

and Qyx ~ ^^'i'Qi'/x- Finally, the set of the three linearly independent third order integrals 
is 



= I -iVo I , = I I , = ( -iPo 






2 

-iP3. 



(6.6) 
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Their linear combination corresponds to the integral L, L = ^{L^"^ + L^^ + L^^). 
The grading operator can be chosen in one of the form^ 

nr \ /7^r o o 

nr , fi = rf = -nr o | . (6.7) 
nrj \ tzTj 

When r is chosen as the grading operator, we have eight nontrivial fermionic integrals of 
motion Sf"^, Sf^, Q]Jxi Qjfxi seven linear independent bosonic integrals S}^ , Qlfx^ 

L^^, L^^ and L^^. In the case of the choice of f as the grading operator, we have nine 
fermionic integrals Sf^ and L"''. The second order integrals Qf^ constitute the set of six 










0\ 









-1 









^0 





I) 





bosonic integrals of motion. Finally, for Fi, we have six bosonic integrals, S"/^, Sf^ and Q^^ 



i/x 

^^volc Cl2 C23 „„A n 

//A' 

and nine fermionic integrals, Sf^, Qjjx^ Qlfx^ We see that the complete set of local 

nontrivial integrals of motion separates into bosonic and fermionic generators in dependence 
on the choice of the grading operator. 

If we start from the first order integrals S*/^, Sf^ and 5/^, their corresponding (anti)- 
commutation relations fl4.1ip and/or (14.1 7p (that depends on the choice of the grading op- 
erator) generate, unlike the N = 2 case, all the six second order integrals Qf/x- '^^^ virtual 
parameter tx here as well as for > 3 can be identified with the shift parameter of one 
of the corresponding subsystems. Each time, however, the intermediate index A in Qf^x 
be changed by employing relation (I3.12p . Anyway, for N >2 the second order integrals Q1/x 
are generated also via the (anti)-commutators of Sf'^ with the third order integrals, see Eqs. 
dOl, (^M, <^M, (O and dEl. 



7 Self-isospectral soliton chains 

Consider now the infinite period limit which produces self-isospectral non-periodic chains. 
It is obtained by putting — )■ 1, when, as we noted, K — )■ oo, 2iK' — )• ivr, and one-gap Lame 
Hamiltonian (12. 2p transforms into that of refiectionless Poschl- Teller system 

In this limit the valence band < -E < A;'^ of the Lame system shrinks into the level 
E = of the unique bound state of the system (17.10 described by the wave function sech x. 
To get a self-isospectral supersymmetric chain of refiectionless Poschl- Teller systems, we 
have to introduce also some restrictions on the displacement parameters. Namely, it is 
necessary to require that all the which appear in the arguments of the chain Hamiltonians, 
Ha = H{x + Ta), should uot go to infinity when A; — )■ 1. By this condition we prohibit that 
in the chain we get in the limit, there could appear free particle systems. 
In such a limit, for x- independent structures we get 

3 (r) coth 2r , Cat coth 2r„fc , gabx coth 2Tab + coth 2Tbx + coth 2Txa , (7.2) 



^There are other possibilities, not reducible to the change of indices and multiplication of the matrix 
elements by —1 (sec the remark at the end of the previous section), which we do not consider here. 
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6ab csch^2rafe , ^abx -(coth 2Ta\ + coth 2rAfe)csch^2rafe . (7.3) 

For the superpotential (12. 8p we have A{x; r) — )■ A^'^{x; r) = tanh(x — r) — tanh(x + r) + 
coth2r. Denoting the hmit of A{xab', Tab) by A^^^, we find 

A^ff = -A^^ = tanh(x + r^) - tanh(x + u) + coth(rb - r^) 

= tanh 2Tab tanh(x + r^) tanh(x + Tb) + 2csch4rab . (7.4) 

For the first order intertwining operator we get 

Vab Xab = ^ ~ ^alf ' ^Ib = -^ba , (7.5) 

and (13. 2p transforms then into 

- XabXba = H^^ + csch22r,b , (7.6) 

where H^^ = H^'^{x + Ta)- The hmit of the second order intertwining operator (13. lip can be 
written in the form that includes in its structure the first order operator = ^— tanh(x+r), 

(P d 

= -AaAl-COt\l2TabXab, (7.7) 

where Aa = Ar^. The first term in the last expression in (17. 7p . unlike the Xab and the x- 
independent multiplier coth 2x^6 in the second term, is well defined when u = Ta ^ Tab = 0. 
This is so because the operator A^-, unlike the limit of the operator D(x;r), is regular for 
r = 0. 

Making use of Eqs. ( 13. 9p and ( 17. 2p . for the limit of the second order intertwining operator 
B we get 

Bab/X -XaxXxb = AaAl + (coth 2r«A + coth 2Txb)Xab = Bab/X , (7.8) 

and find that the limit of the third order integral can be presented in terms of the introduced 
first order operator A^-, 

Va^-Za, where = . (7.9) 

Since the first order operator Xab as well as the second order operators (17. 7p and (17. Sp 
intertwine the Hamiltonians H^^ and H^^, the operator |30j 

Yab = AaAl (7.10) 

is also the intertwining operator, YabH^^ = H^^Yab- Operator (I7.10p corresponds to the 
infinite limit of the virtual displacement parameter, — )■ oo (or, tx — > — cxd), applied to 
Bab/X, see Fig. 4. 

The intertwiner (I7.10p . unlike (17. 7p and (17. Sp . is regular for Tab = (r^ = t^), when it 
reduces just to H^^ , 

Yaa = AaAl = ■ (7.11) 
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Another product of the same operators produces (for any value of the parameter Ta) the free 
particle Hamiltonian shifted for an additive constant, 

4A = -^ + l = ^o. (7.12) 

In accordance with (17. lip and fl7.12p . the first order operators Aa and intertwine the 
Poschl- Teller system with a free particle, AlH^^ = HqAI, H^^ Aa = AaHo. From (17.41) we 
find that if r;, — oo while is kept to be finite, A^ff reduces to tanh(x + ra). In such a limit, 
H^^ reduces to Hq, = Hq, Xab transforms into Aa, while the second order operators 
(17. 7p . (17. 8p and (I7.10p transform into Aa-^ and linear combinations of this operator and the 
first order operator Aa- The second order operator we have gotten intertwines Hj^^ with 
Hq, Aa-^Ho = H^^Aaj-, but this relation produces nothing new since it is a consequence 
of the conservation of ^ for a free particle system Hq, -^Hq = Hq-^, and of the already 
known relation H^^Aa = AaHo. On the other hand, the product of this operator with 
the intertwiner Al, that acts in another direction between Hq and H^^, shows that the 
Lax operator of the Poschl- Teller system, Za, is nothing else as the Darboux-dressed free 
particle momentum [51]. From (17. 8p one can get a relation which involves the first order 
intertwiners Aa and Xab, and the free particle momentum. Taking the limit r;, — )■ oo in 
both representations for Bat/x, we get the relation Aa — coth2raA) = XaxAx , and also 
its conjugate, A{Xxa = {-£ + coth2raA) A^. 




Figure 4: The first order Darboux displacement generator Xab translates the eigenfunctions 
of H^^ into those of H^^. The second order intertwiner Yab makes the same via the virtual 
free particle system that is a Poschl- Teller system translated to infinity. Lax integral Za is 
presented as a Darboux dressed form (17.90 of the free particle integral d/dx. 



The described limit apphed to the chain (14. ip gives a non-periodic self-isospectral chain 
described by the Hamiltonian 

H^^ = diag{H^^,...,H^'^). (7.13) 

The nontrivial integrals of such a system are given by the infinite period limit applied to 
the integrals S'f^, ^'^d t'Y employing relations (17. Sp . (17. 8p and (17. 9p . As in the 
periodic case, the set of the corresponding third order integrals Lf^""^ contains only in- 
dependent integrals. Instead of QfK"'^, one can work with the second order matrix integrals 
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Qf = limT-;^^oo Qi/x \ which are obtained by the change of Bab/x for Yab- The corre- 
sponding superalgebra generated by these integrals, as in the periodic case, depends on the 
choice of the grading operator, and its concrete form can be computed by making use of the 
relations presented above. Again, we get a closed non-linear superalgebra, the nonlinearity 
of which originates from the polynomial dependence of the superalgebraic structure func- 
tions on the Hamiltonian "H^^, which plays a role of the multiplicative central charge. Since 
we have defined Qf^"^ as the Qf/J^"'^ taken with the same virtual parameter ta = oo, the 
set of the second order integrals Qf"^""^ together with the Hamiltonian form the closed 
nonlinear sub-superalgebra. 

If some of the displacement parameters are taken to be infinite, corresponding Hamil- 
tonians H^^ transform into those of the free particle. In this case we loose the property of 
the self-isospectrality of the chain, but the supersymmetric structure is still present and can 
also be computed by employing the relations discussed above. 



8 Discussion and outlook 

We investigated the unusual nonlinear supersymmetric structure of the self-isospectral crys- 
talline chains formed by the arbitrary number > 2 of mutually displaced periodic one-gap 
Lame systems, and of the associated non-periodic self-isospectral soliton chains described 
by reflectionless Poschl- Teller Hamiltonians. It is generated by the A^(A^ — 1) integrals of 
motion which are the first order differential matrix operators, by the same number of the 
second order matrix integrals, and by the third order Lax integrals. The supersymmetry 
admits distinct choices for the grading operator, that classifies these integrals as bosonic and 
fermionic operators in different ways. For instance, for the example of the A^ = 3 chain from 
Section 6, three choices of the grading operator indicated in Eq. (16. 7p classify the complete 
set of 15 nontrivial local integrals of motion as, respectively, 7 + 8 = (2 + 2 + 3) + (4 + 4 + 0), 
6 + 9 = (0 + 6 + 0) + (6 + + 3) and 6 + 9 = (4 + 2 + 0) + (2 + 4 + 3) bosonic+fermionic 
generators, where the first, second and third numbers in each parentheses correspond to the 
numbers of the first, second and third order differential matrix operators. In dependence on 
the chosen grading operator, one of the third order integrals, L, is the bosonic central charge, 
or the fermionic supercharge to be a square root of the matrix spectral polynomial of the 
A^-chain In the latter, unlike the former, case the spectral polynomial appears explicitly 
just in the superalgebraic relations. This reveals the identifying characteristic of the Lax 
integrals: they recognize the band edge states in the spectrum of each Lame subsystem by 
annihilating them. Another peculiarity is that the set of all the second order integrals of 
motion taken with the same virtual parameter generates together with the Hamiltonian a 
nonlinear sub-superalgebra. 

The lowest A^-fold degenerate energy level was chosen to be zero, and the spectra of the 
self-isospectral chains described by the second order matrix Hamiltonians do not depend on 
the values of the displacement parameters. On the other hand, according to Eqs. (I4.10p and 
(12.1 op . the spectra of the first order matrix integrals Sf^ depend on the mutual shifts Tab, and 

•^As we noted at the very beginning, this happens even in the case = 1 for the unextended Lame 
system (|2.2p . which is described by a hidden, bosonized supersymmetry with the reflection TZ identified as 
the grading operator [T^l H] . 
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blow up when Tab tends to zero (modulo the period in the case of the crystalline chain) . This 
indicates on another possibility to interpret the systems by identifying a suitable combination 
of the first order integrals as a Hamiltonian. For instance, for N = 2, one can treat 5*^ as the 
integral T-L{i), or, for N = 2n we can choose 7{(i) = 5*^ + Sf^ + . . . ^^"-isn^ r^^^ lower index 
indicates that the Hamiltonian is of the first order (Dirac) nature, which can be considered 
as a kind of Bogoliubov-de Gennes Hamiltonian in the Andreev approximation. Such "H^i) in 
the N = 2 case was considered, for instance, in the physics of conducting polymers [T5ll35] . or 
as a Hamiltonian that describes the kink-antikink crystal |25l [26] (or, kink-antikink baryons 
in the non-periodic limit case [36]) in the Gross-Neveu model. Therefore, the N > 2 chains in 
such a reinterpretation with the first order Hamiltonian would provide some generalization 
of the known N = 2 models, in which spectral gaps are governed by the displacement 
parameters of the corresponding second order chains. The interesting peculiarity of such first 
order systems is that they possess the own nonlinear super symmetry. Indeed, the operator 
Fi = TZTT, where F is given by Eq. fl4.9p . commutes with T-L(i), and can be identified as 
a grading operator for such a first order system. The Lax operator L anti-commutes with 
Fi, and the latter classifies 'H(i) and L as, respectively, bosonic and fermionic generators. 
Since L commutes with it can be considered as a fermionic supercharge, whose square, 
in accordance with Eqs. fl4.10p and fl5.6p . gives some polynomial of order six in 'H(i) q For 
n > 1, the system 'H(i) has also other nontrivial integrals of motion, see Eq. 04.240 with 
/ = 1. Such a nonlinear supersymmetry in the first order system T-L{i) could not be revealed, 
however, with the choice of the grading operator in the reflection-independent form (14.90 
which identifies the Lax integral L as the bosonic operator and ^-(i) as the fermionic one. 

Another interesting possibility for generalization of the results is to identify some linear 
combination of the second order matrix operators as a Hamiltonian, for instance, by taking 
H(2) = Q\% + Ql% + - ■ ■ Ql'Jx^^'' in the case of = 2n. The spectrum of H(2) like that of 
depends on the shift parameters, see (I4.12p . For A^ = 2 or A^ = 4, such a second order matrix 
Hamiltonian has a nature to be similar to that of the Hamiltonian in the physics of bilayer 
graphene [2Z] ■ With respect to the grading operator F = TZT, the Hamiltonian 7/(2) and the 
third order operator L are, respectively, the bosonic and fermionic operators. The operator 
L as well as the operators L\^, . . ., Ll^~^'^^, see Eq. (14.260 with / = 1, are the supercharges 
of the nonlinear supersymmetry of the system described by the unusual second order matrix 
Hamiltonian 'H(2). Notice a special role played by the reflection-dependent grading operator 
F = TZT for revealing the supersymmetric structure in the indicated unusual second order 
system 'H(2)- 

The peculiarity of the non-periodic case in comparison with the periodic one is that the 
chain subsystems can be related there, particularly, by the second order intertwiners (l7.1Up . 
As the intermediate (virtual) system in this case, there appears a free particle system. The 
latter can be treated as the Poschl- Teller system (17.10 . = H['^{x + ta), displaced to 
infinity, ta — t- oo. It is due to such a relation the Lax operator (17. 9p has a nature of a dressed 
free particle momentum operator, and eigenstates of H^^ can be obtained by the Darboux 
transformation of the corresponding free particle eigenstates. We have, unfortunately, no 
such simple relation with a free particle in the periodic case. 

We considered the case of self-isospectral Hermitian chains with real displacements. The 

Such a nonlinear supersymmetry in the first order systems 'H(i) was discussed in [221130) for the simplest 
case of = 2 chains; it appears particularly in the twisting of carbon nanotubes, see [38j . 
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construction can be generalized for the case of complex shift parameters. The correspond- 
ing supersymmetric structure can be interesting then in the context of the physics of PT- 
symmetric systems [221110], where, again, the discrete transformation operators, particularly 
spatial reflection, prove to play a fundamental role [¥T] . 
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